Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 

Updated Measurement of B{B, Di*^^Di*^~) and Determination of Ar, 
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Using fully reconstructed Bs mesons, we measure exclusive branching fractions for the decays Bs —>■ -Ds*'^£'s*' • 
The results are B{B0^D+D-) = (0.58l[;;J^ ±0.13)%, B{BO^D*^Df) = (1.8 ±0.2 ±0.4)%, and B(B°^ 

Dl+Dt-) = (2.0 ±0.3 ±0.5)%; the sum is B{B° ^ oi*'''^ d[''>~ ) = (4.3 ±0.4 ± 1.0)%. Assuming these decay 
modes saturate decays to CP-even final states, the branching fraction determines the relative width difference 
between the Bs CP-odd and CP-even eigenstates. Taking CP violation to be negligibly small, we obtain 
Ar^/r^ = 0.090 ± 0.009 (stat.) ± 0.022 (syst.), where is the mean decay width. The results are based on a 
data sample collected with the Belle detector at the KEKB e+e~ collider running at the T(55) resonance with 
an integrated luminosity of 121.4 fb~^. 



,3 1. Introduction 

^ An e^e collider running at a center-of-mass (CM) energy corresponding to the T(5S') resonance can produce 

significant amounts of Bs*''^!*'' pairs [T| H]. The Belle detector [3] at the KEKB asymmetric-energy e+e~ 
collider [4] has collected a data sample corresponding to an integrated luminosity of 121.4 fb~^ at the T{bS) 
resonance {s/s = 10.87 GeV). This sample has allowed us to make the most precise measurement of the branching 

^""^ fractions of i3° —> decays [S]. These Cabibbo-favored final states are expected to be predominantly 

CP-even [5] and dominate the difference in decay widths AF^^ between the two B^-B^ CP eigenstates [B]. 

^ We report preliminary results of the updated branching fraction measurements, which replace our previous 

D measurement based on 23.6 fb~^ of data H. 
I 

^ 2. Event Selection 

I The Belle detector consists of a silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), 

J> an array of aerogel threshold Cherenkov counters (ACC), time-of- flight scintillation counters (TOF), and an 
electromagnetic calorimeter (ECL) comprising of CsI(Tl) crystals located inside a superconducting solenoid coil 
that provides a 1.5 T magnetic field. For charged hadron identification, a likelihood ratio is formed based on 
a dE/dx measurement in the CDC and the response of the ACC and TOF. Only good quality charged tracks 
originating from near the e+e" interaction point are accepted. For charged kaon tracks, this likelihood ratio 
is required to be > 0.60; the tracks not satisfying this requirement are identified as pions. The kaon likelihood 
requirement is ^90% efficient and has a tt* misidentification rate of ^10%. With the exception of the tracks 
I originating from decays, low-momentum charged tracks with P < 100 MeV/c are rejected. Neutral Kg 
• • candidates are reconstructed from tt^tt" pairs having an invariant mass within 10 MeV/c^ of the nominal Kg 
^ ^ mass [S] and satisfying momentum-dependent requirements based on the decay vertex position [Sj . 
K> Neutral tt^ candidates are reconstructed from 77 pairs having an invariant mass within 15 MeV/c^ of the ir'^ 

mass with each photon having a laboratory energy greater than 100 MeV. Neutral and charged K* candidates 
are reconstructed from a K and 7r+ having an invariant mass within 50 MeV/c^ of Mj^,. Neutral candidates 
are reconstructed from K~^K^ pairs having an invariant mass within 12 MeV/c^ of M^. Charged p'^ candidates 
are reconstructed from tt+tt^ pairs having an invariant mass within 100 MeV/c^ of M^+. 

We reconstruct Df candidates using six final states: 4>tt+ , K^K+, K*°K+, (t>p+ , K^ K*+ , and K*^K*+. 
The invariant mass windows used are 10 MeV/c^ (~ 3cr) for the three final states containing K* candidates, 
20 MeV/c^ (2.8cr) for (t>p'^ , and 15 MeV/c^ (^ Aa) for the remaining two modes. For the three vector- 
pseudoscalar final states, we require |cos0jj(,[| > 0.20, where the hehcity angle 9^^^ is the angle between the 
momentum of the charged daughter of the vector particle and the direction opposite the Dg momentum, in the 
rest frame of the vector particle. 

We combine candidates with photon candidates to reconstruct -D*^ — > -D^7 decays, and we require 
that the mass difference ^f^+^ ^ ^^nt within 12.0 MeV/c^ of the nominal value (143.8 MeV/c^), where 

Df denotes the reconstructed candidate. This requirement (and that for the Df mass) is determined by 
optimizing a figure-of-merit S/ \/S + B, where S is the expected signal based on Monte Carlo (MC) simulation 
and B is the expected background as estimated from £)+ sideband data. We require that the photon energy in 



2 
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the CM system be greater than 50 MeV, and that the energy deposited in the central 3x3 array of cells of the 
ECL cluster contain at least 85% of the energy deposited in the central 5x5 array of cells. 

We select B° Dl+Dl^ , Dl^Df, and DfD~ decays using two quantities evaluated in the e+e^ CM 
frame: the beam-energy-constrained mass My^^ = \/ ^bcam ~ P% ' ^"^"^ energy difference AE — Eg — i^beam' 
where and E^ are the reconstructed momentum and energy of the candidate, and -Ebcam is the beam 
energy. We determine our signal yields by fitting events in the region 5.25 GeV/c'^ < M^^^ < 5.45 GeV/c^ 
and -0.15 GeV < AE < 0.10 GeV. Within this region, the modes T{5S)^BJ3^, BJ3* and B*B* are well- 
separated as the 7 from B* -^B^^ decay is not reconstructed. We expect only small amounts of signal in B^B^ 
and BgB* and thus do not use these modes for the branching fraction measurement. In order to simplify the 
fit, we fix their relative ratios, which are determined from fully reconstructed — >■ D~Tr~^ decays as described 
in Ref. [H]. 

When multiple -6° — > D^*''^ D^*"^ candidates are reconstructed in an event, we select the candidate that 
minimizes the quantity 
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where AM — Mjj, — Mj^ , the quantities M^ and AM are reconstructed, and the summations run over the 
two daughters and the possible Z3*+ daughters {N = 0,1,2) of a _B" candidate. The mean mass Af^ 
and widths a^^ and cr^j^j are obtained from MC simulation and calibrated for any data-MC difference using a 

B'^ ^ oi*^^ D~ sample in 563 fb~^ of data at the T{4S) energy. Approximately half of the events have multiple 
candidates according to MC simulation, and this criterion selects the correct candidate 83%, 73%, and 69% 
of the time for DfD~, Dl'^Df, and -D*+Z)*^ final states, respectively. 

The background from e'^e~ ^qq {q = u, d, s, c) continuum events is rejected using a Fisher discriminant based 
on a set of modified Fox- Wolfram moments [TO]. This discriminant distinguishes jet-like qq events from more 
spherical B^^s^B^^-^ events, and is used to calculate a likelihood £s i^gq) for an event assuming the event is signal 
{qq background). We require the ratio TZ = Cs/{Cs -\- Cqq) to be >0.20. This selection is 93% efficient for signal 
and removes > 62% of qq background. The majority of the background consists of T [bS) ^ B^f^ B^a^ ^ D'^ X; 
T{5S)^BBX (where bb hadronizes into 5°, 5°, or B±); and B^-^ D^j {2317) Di*\ B,^ Dfj{2m)Di*\ and 
B^-i' DfDfir" decays. The last three processes peak at negative values of AE, and their yields are expected 
to be very small assuming their branching fractions are similar to analogous B^^D^jD^*^ decays. 

Signal yields are measured using a two-dimensional extended unbinned maximum-likelihood fit to the M^^-AE 
distributions. For each signal decay, we include probability density functions (PDFs) for signal and background. 

We use a single PDF for background which consists of qq, B^*^ Bi*^ — >• D'^X, and T(5S') -> BBX events. 
The background PDF is constructed using an ARGUS function [TT] for Mj^^, and a second-order Chebyshev 
polynomial for AE. The two parameters of the Chebyshev function are taken from the data in which one of the 
Ds "candidates" is required to be within the mass sideband. The signal PDFs have three components: correctly 
reconstructed (CR) decays; "wrong combination" (WC) decays in which a non-signal track or photon is included 
in place of a true daughter track or photon; and "cross-feed" (CF) decays in which a Dl^Df or Dl^D*~ is 
reconstructed as a DfDj or D^^DJ , respectively, or else a DfD~ or Dl'^Df is reconstructed as a Dl'^Df or 
D*'^D*~ . For these CF candidates AE is shifted by 100-150 MeV, but M^^ remains almost unchanged. When 
the 5° is not fully reconstructed, e.g. due to losing the 7 from _D*+— >-I?+7 (CF-down), a negative shift in AE 
is observed. Conversely, in the case where the signal decay has gained a photon (CF-up), AE is typically shifted 
higher. The PDF for CR events is modeled with a Gaussian for M^^ and a double Gaussian with common mean 
for AE. CF and WC events have more complicated distributions. All signal shape parameters are taken from 
MC and calibrated using B'^^-^oi*^ n+ and B"^Di*^^D- decays. The fractions of WC and CF-down events 
are taken from MC simulation. The fractions of CF-up events are difficult to calculate accurately from MC 
simulation as not all partial widths are measured; thus they are allowed to vary in the fit. As the CF-down 
fractions are fixed, the three distributions {DfD^, D'*^Df, and -D*+-D*^) are fitted simultaneously [12]. The 
CF fractions are listed in Table H] 



2.1. Results 

We measure t' 

fraction fg,^, — (87.0 ± 1.7)% [TO] . The fit results are listed in Tablepl] and projections of the fit are shown 



We measure the signal yields for B'^ oi*^^ oi*^ decays using 7.1 ± 1.3 million Bi*^ B^*'' pairs with a B*Bl 
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Table I: Fractional distribution of the signal reconstruction types from MC simulation of B° decay modes. 



B° Mode 


RC 


WC 


CF I 




CF II 


d+d: 


76.1 


6.0 fixed 


17.1 {-^DI^DJ) 


0.8 ( 






44.4 


38.5 fixed 


8.2 (-^D+Z)7)fixed 


8.9 ( 


-^D*+D*-) 


Dl+Dl~ 


31.8 


37.6 fixed 


2.0 {~>DtD-) fixed 


28.6 (- 


>D*^Df) fixed 



in Fig. g 



The branching fraction for channel i is calculated as B.^ = Y^/(e' 
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2), where 



is the fitted CR yield, and s\j(j is the MC efficiency with intermediate branching fractions |H] included. 
The efficiencies £\f(j include small correction factors to account for differences between MC and data for kaon 
identification. 

The systematic errors are listed in Table III The error due to PDF shapes is evaluated by varying shape 
parameters by zblcr. The systematic error for the fixed WC and CF-down fractions is evaluated by repeating 
the fit with each fixed fraction varied by ±20%. The uncertainties due to identification and tracking are 
^ l%(momentum-dependent) and 0.35% per track respectively. As the longitudinal polarization fraction (/^) 
of — >• Dl'^Dl~ is not measured yet, we assume to be the world average (WA) value for the analogous 
spectator decay D'*^D*^: 0.52 ± 0.05 [8]. The related systematic error is taken as the change in B when 
is varied by twice the error on the WA value. Significant uncertainties arise from branching fractions and 
the fraction of T(5S') decays producing B^ mesons, which are external factors that are expected to be measured 
more precisely in the future. The statistical significance given in Table [ll| is calculated as 2 ln(£o/'Cmax), 
where Cq and £max are the values of the likelihood function when the signal yield is fixed to zero and when 
it is set to the fitted value, respectively. We include systematic uncertainty in the significance by smearing the 
likelihood function by a Gaussian having a width equal to the total systematic error obtained for the signal 
yield. 
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Figure 1: Af^,, and AE projections of the fit result. The columns correspond to B'^ — )■ D'^ DJ (right), B^' — >■ Dl^Df 
(middle), and B^ D^^ Dg~ (left). The red dashed curves show RC-fWC signal, the blue-purple solid curves show CF, 
the grey solid curve shows background, and the black solid curves show the total. 
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Table II: Signal yield {Y), efficiency including intermediate branching fractions (e), branching fraction (B), and signal 
significance (S) including systematic uncertainty. The first errors Hsted are statistical, and the second are systematic. 



Mode 


Y 


e 


B 


S 




(events) 


(xlO-*) 


(%) 




DtD- 


33.lt|« 


4.72 


0.58lg;Jg ±0.13 


11.5 




U.5tli 


2.08 


1.8 ± 0.2 ±0.4 


10.1 


did: 




1.01 


2.0 ±0.3 ±0.5 


7.8 


Sum 






4.3 ± 0.4 ±1.0 





Table III: Systematic errors (%). Those listed in the top section affect the signal yield and thus the signal significance. 



Source 


DfDj 






D^D^ 




Dt+Df 






±cr 


— CT 


±cr 




— O" 


±cr 


— O" 


Signal PDF Shape 


2.7 


2.2 


2.2 




2.4 


5.1 


3.8 


Background PDF Shape 


1.5 


1.2 


1.3 




1.4 


2.9 


2.2 


WC ± CF fraction 


0.7 


0.6 


4.6 




4.5 


6.2 


6.2 


TZ requirement {qq suppression) 


3.1 


0.0 


0.0 




2.7 


0.0 


2.1 


Best candidate selection 


5.5 


0.0 


1.5 




0.0 


1.5 


0.0 


Identification 


7.0 


7.0 


7.0 




7.0 


7.0 


7.0 


Ks Reconstruction 


1.1 


1.1 


1.1 




1.1 


1.1 


1.1 


7r° Reconstruction 


1.1 


1.1 


1.1 




1.1 


1.1 


1.1 


7 






3.8 




3.8 


7.6 


7.6 


Tracking 


2.2 


2.2 


2.2 




2.2 


2.2 


2.2 


Polarization 


0.1 


0.1 


0.8 




0.7 


0.5 


1.0 


MC statistics for e 


0.2 


0.2 


0.4 




0.4 


0.5 


0.5 


_Di*' Branching Fractions 


8.6 


8.6 


8.6 




8.6 


8.7 


8.7 










18.3 








f SJB* 








2.0 








Total 


22.7 


21.8 


22.6 




22.8 


24.6 


24.3 



In the heavy quark limit vifith (m^ — 2m^) — >■ and — >■ oo, the dominant contribution to the decay width 

comes from oi*^^ ui*^ decays |6l[l4j. Assuming negligible CP violation, the branching fraction is related 
to Ar^ as Arjr, ^2B/{1-B). Inserting the total B from Table |!T] gives 



AT. 



0.090 ± 0.009 ± 0.022, 



(2) 



where the first error is statistical and the second is systematic. This result is in good agreement with the current 
WA |S] and is consistent with theory [T5]. There is a theoretical uncertainty arising mainly from the CP-odd 
component in B'^ D*~^ D*~ and the unknown contribution of 3-body final states. 
If a CP-violating phase (j)s is allowed, the above relation becomes 



4B{B'^^Di*^Di*y)^(^) 

\ cos (f)s ) 



1 



COS ( 



1 — COS ( 



l±Ar?^ 1-Ar?^ 



(3) 



where (/-^ = Arg{M^2l^\2) [IS]- Fig. [2] plots AF^ as a function of 0^ for our measurement. 
In summary, we have measured the branching fractions for P^ — )■ £'i*^"''Z)i*^ using e^ e~ data taken at the 
T(5S') resonance. Our results constitute the first observation of P° — )• D*'^D*^ (Scr significance). Using the 
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Figure 2: The width difference AF^ as a function of (ps- One-sigma band and SM value are shown for comparison. 



total measured branching fraction 6(5° ^ Z^i*' ) = (4.3 ± 0.4 ± 1.0)% and assuming no CP violation, 
we determine the relative BJ3^ decay width difference to be AF^/F^ = 0.090 ± 0.009 ± 0.022. 
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